Priedel oscillations, impurity scattering and temperature dependence of resistivity in 
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We show that Friedel oscillations (FO) in grapehene are strongly affected by the chirality of elec- 
trons in this material. In particular, the FO of the charge density around an impurity show a 
faster, 8p ~ r -3 , decay than in conventional 2D electron systems and do not contribute to a linear 
temperature-dependent correction to the resistivity. In contrast, the FO of the exchange field which 
surrounds atomically sharp defects breaking the hexagonal symmetry of the honeycomb lattice lead 
to a negative linear T-dependence of the resistivity. 

PACS numbers: 73.63.Bd, 71.70.Di, 73.43.Cd, 81.05.Uw 
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Screening strongly influences properties of impurities 
in metals and semiconductors. While Thomas- Fermi 
screening suppresses the long-range tail of a charged im- 
purity potential, Friedel oscillations (FO) of the electron 
density around a defect [l[ are felt by scattered electrons 
at a distance much longer than the Thomas- Fermi screen- 
ing length. Friedel oscillations originate from the singular 
behavior of the response function of the Fermi liquid at 
wave vector 2kp. At zero temperature the decay of the 
amplitude Sp of these oscillations with distance r from 
the impurity obeys a power law dependence. In a non- 
relativistic degenerate two-dimensional (2D) Fermi gas 
2], bp oc cos(2kpr + S) / r 2 . In 2D electron systems Bragg 
scattering off the potential created by these long-range 
FO strongly renormalises the momentum relaxation rate, 
r _1 for quasi-particles near the Fermi level, e ~ ep, 
which leads to a linear temperature dependence of the 
resistivity 0, 0, H, @] m a 'ballistic' temperature range 
ep > T > h/r confirmed in recent experiments on semi- 
conductor heterostructures and Si field-effect transistors 

0- 

Graphene-based transistor || Q is a recent invention 
which attracts a lot of attention. Improvement of the 
performance of this device requires identification of the 
dominant source of electron scattering limiting its mo- 
bility. Those can be structural defects of graphene lat- 
tice (vacancies and dislocations), substitutional disorder, 
chemical deposits on graphene and, importantly, charges 
trapped in or on the surface of the underlying substrate. 

In this Letter we investigate the effects of screening of 
scatterers in graphene, and show how the phenomenon of 
FO can be used to gain insight into the microscopic na- 
ture of disorder. Graphene (a monolayer of graphite) is a 
gapless 2D semiconductor with a Dirac-like dispersion of 



carriers 



HI El 



In this material, it is the Thomas- Fermi 



screening which is responsible [12[ for the experimentally 
observed linear dependence of graphene conductivity on 
the carrier density 0, 0|. Moreover, quasiparticles in 
graphene possess chiral properties related to the sublat- 
tice composition of the electron wave on a 2D honeycomb 



lattice Below we show that, due to the latter pe- 

culiarity of graphene, FO of the electron density decay 
as Sp ~ r~ 3 at a long distance from an impurity [faster 
than in a usual 2D metal] and the linear temperature- 
dependent correction to the resistivity, 



R(T) - R(0) = 



(1) 



is caused only by the FO of the exchange field and is 
strongly sensitive to the microscopic origin of scatters. 
In Eq. ([1]), r" 1 is the backscattering rate specifically 
from atomically sharp defects distorting the hexagonal 
symmetry of the honeycomb lattice: structural defects, 
chemical deposits and substitutional disorder. Note that 
Coulomb scatterers do not contribute towards this result. 
Thus, we predict that in graphene-based devices where 
scattering is dominated by charges trapped in the sub- 
strate (or on its surface) the temperature dependent re- 
sitivity, [R(T) - R(0)]/R ~ (T/e F )(r/r») is much weaker 
than in conventional 2D semiconductor structures. We 
propose to use the measurement of the resistivity in a 
temperature range £f > T > ft/r as a probe of micro- 
scopic composition of disorder. 

Electrons in graphene can be described using 4- 
component Bloch functions [4>k + ,A> 4>~k + ,b, 0k_ b, 
<Pk-,a]i which characterise the electronic amplitudes on 
the two sublattices (A and B) and valleys K + and K . 
Since for the "ballistic" temperature regime Tr/h > 1 in- 
terference between waves scattered from FO around dif- 
ferent impurities can be neglected, we study screening of 
an individual impurity. The single-particle Hamiltonian 
in the presence of one impurity is 



-ifeE-V + u5(r). 



(2) 



In this expression X = (Y* x ,Y, y ) is a two-dimensional 
vector whose components are 4x4 matrices E x = H z ®cr x 
and £„ = IL Z (8 <J y . Together with E 2 = IIo ® <j z they 
the generator algebra of the unitary group SUf . 
and cT x ^y yZ are sets of three Pauli matrices 



form [l3l 
Here II 
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acting on the valley and sublattice indices, respectively, 
and Ilo (<7o) are 2x2 unit matrices acting on the valley 
(sublattice) spaces. Another set of three matrices, A x — 
Ihc <8> cr z , A y = Hy <x> a z and A z = U z <8 er , such that 
[A h ,A h ] = 2ie hhl A h satisfies [£ S ,A Z ] = 0. Therefore, 
they commute with the first term in H sp and generate 

the unitary group SUj = {e lbnA } describing the valley 
symmetry of the Dirac Hamiltonian in graphene. 

For a non-magnetic defect the matrix u in Eq. Q 
should be hermitian and time-reversal symmetric. It can 
be par ametrised using 10 independent real parameters 
13, 111. One can check [3, EH that all the operators T, s 
and A; change sign under the time-reversal transforma- 
tion, so that the 9 products S S A; are t — > —t invariant 
and together with the unit matrix I = IIo(8><7o can be used 
as a basis to represent a non-magnetic static disorder: 

u = ul + u si^sAi =wi + A z £ z + X-£, (3) 

s,l—x,y,z 

X s = ^ u siAi. 

The term ul in Eq. ([3]) represents an electrostatic poten- 
tial averaged over the unit cell. We attribute this term 
to charged impurities with a sheet density n c . Various 
atomically sharp defects [l6| with the 2D density ridei 
can be parametrized using nine independent real parame- 
ters u s i and time- inversion-symmetric [la ] matrices £ s Aj. 
In particular, u zz describes different on-site energies on 
the A and B sublattices. Terms with u xz and u yz take 
into account fluctuations of A ^ B hopping, whereas 
u sx and u sy (s = x, y, z) generate inter-valley scattering 
(whose presence has been revealed tl3[ by the observation 
of weak localisation in graphene [171]). 

Thomas-Fermi screening. The momentum relaxation 
of chiral electrons in graphene due to scattering off 
Coulomb impurities is determined by the anisotropic dif- 
ferential cross-section, 

^)~<si„ (e /2) cos 2 (9/2). 

The factor cos 2 (0/2), where 9 is the scattering angle, re- 
flects the absence of back-scattering of graphene elec- 
trons off the electrostatic potential 



density n e [lz 
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In the 

Thomas-Fermi approximation, the Fourier transform u q 
of the electrostatic potential of a single charged impurity 
screened by 2D electrons [1(3] is u q = 2ir(e 2 /x)/(q + K )> 
where k = 8ir'je 2 /x = ^kpr s ; 7 = kp / (2whv) is the den- 
sity of states per spin and valley, x is the dielectric con- 
stant, and r s = e 2 /x?iv (for a graphene sheet placed on 
Si substrate, r s ~ 1). 

In a structure with sheet density n c of charged im- 
purities, the resistivity, R — [2e 2 jv 2 r]~ 1 is determined 
by the momentum relaxation rate, r _1 ~ n^^hr 1 ((1 — 
cos9)w(9))g, and it is inverse proportional to the carrier 



R 



h n c Airri(r s )r 2 
4e 2 n e (1 + 4r s ) 2 



(4) 



where r/(r s ) is a monotonous function such that 77(0) = 1, 
77(1) « 0.3 and 17(00) = j. It has been established that 
a contribution from atomically sharp disorder towards 
R depends on the Fermi energy only logarithmically 14 1 . 
Therefore, the empirical relation n e R — const established 
in the experimentally studied graphene structures with 
high carrier densities [1, @| indicates that charges located 
in the underlying substrate or on its surfaces are the main 
source of scattering. In Eq. ((3]) it is represented by ul. 

Analysis of Friedel oscillations. Below, we analyse the 
FO of the density matrix of 2D electrons in graphene sur- 
rounding various types of scatterers, Eq. ([3|). The plane 
wave eigenstates of the Dirac Hamiltonian — ihvll-V are 
V'k^(r) = e lkr |k£). Here, the spinor |k£) has a definite 
valley projection, A 2 |k^) = £|k£), £ = ±1 and chirality: 
the projection £ • n = ±1 of the operator £ on the di- 
rection of motion, n = k/fc. Below, we describe a chiral 
electron with £ • n = 1 using the polarization matrix 



«=±i 



|k£)(k£| = ±(l + E-n), 



(5) 



Due to scattering off the impurity, the plane wave 
eigenstates of the Hamiltonian acquire a correction Sipk^ , 
which in the Born approximation is given by 



^k,j(r) 



d 2 p 



(2tt) 2 E(k) - hvY, ■ p + iO 



u|k£). (6) 



Here, E(k) — ±hvk, where the sign +(— ) corresponds 
to the conduction(valence) band states. For graphene 
with the Fermi energy ep positioned in the conduction 
band, and for T < ep, we shall disregard the valence 
band contribution. Then, the correction to the electron 
density matrix induced by the impurity is 

5p(r,r') = /(^^( fc WMW®^(r')+^. (7) 

The matrix 6p(r, r') contains a slowly decaying oscil- 
latory part, which is due to a jump in the Fermi function 
np(k). After substituting Eq. © into Eq. ([7]) for the den- 
sity matrix we find that in the region 1 < kpr < ep/T, 
the oscillating part of the 'local' density matrix, p(r, r) 
calculated to leading order in 1/r is given by 



dp(r,r) 



kp 



8n 2 v 2ir 2 



■ s n us — n 



h.c. 



(8) 



s n ws_ n = \(X Z + in x X)(£ a — in X £), 



where n = r/r and a x b = a x b y — a y b x . The matrix 
structure of Sp describes the distribution of charge at the 
atomic length scale: oscillations with the wave vector 
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K = K + — K_ (related to h. XtV ) superimposed with the 
oscillations between the A and B sublattices (related to, 
e.g., E a ). 

The FO in the density matrix © do not lead [22J to 
oscillations in the charge density, Sn e (r) = Tn5/5(r,r) be- 
cause for any pair of matrices S s and A;, Tr£ s A; = 0. To 
find Sn e we evaluated Sp up to order r~ 3 in the l/(k F r) 
expansion and arrived at the leading non-vanishing con- 
tribution to the FO of the electron density coming from 
the diagonal disorder, ul 



5n e (r) un e cos(2fcpr) 
n e ep (2fcpr) 3 



k F r > 1, 



(9) 



where for a Thomas-Fermi screened impurity charge 
(un e /e F ) ~ 2r s /(l + 4r s ) - 1. The FO in Eq. © de- 
cay with the distance from the Coulomb impurity faster 
than in a conventional 2D electron gas (where FO obey 
the 1/r 2 law) . This is due to the absence of backscatter- 
ing of a chiral electron off the potential ul conserving the 
sublattice state [Hj]. As a result, Bragg scattering off FO 
formed around the Coulomb impurity is suppressed and 
does not lead to a linear T-dependence of the resistivity. 

Interaction correction to resistivity. Since atomically 
sharp defects do not generate FO in the electron density, 
they also do not induce an oscillating Hartree potential, 
though they generate a non-local exchange field. Below 
we investigate using the Born approximation how the ex- 
change field created by the FO in Eq. renormalises 



impurity scattering and leads to a linear T-dependence of 
the resistivity. We write down the Hartree- Fock potential 
created by the impurity as a sum Hn + Hp, 

H H = 2 J dvdv'V(v - r')Trp(r',r / )V' t (r)i/'(r), 
Hf = - J drdr'V{r -r')^(r)p(r,r')iP(r'), 

where V(r) is the (screened) Coulomb potential averaged 
over the unit cell [21]. Here, we suppressed the electron 
spin indices, for brevity. The factor of 2 in Ha is the 
result of summation over spin channels. 

The leading correction Sw(p, 9) to the differential cross- 
section w(6) of scatteringoff the screened defect is a re- 
sult of the interference [3j between the wave scattered 
by the defect itself and the wave Bragg-reflected by the 
Hartree, Fig. [TJa) and the Fock, Fig. HJb) potentials of 
the FO, The resulting correction to the momentum re- 
laxation rate is 



d9 



n dcf v I dpn' F (p) I — (1 - cos6)6w(p, 9) 



After substituting the Born amplitude, A(9) of the scat- 
tering process of the defect itself and the amplitudes, 
6A a ( b )(0) of the two processes in Fig. 1 (a) and (b) into 
the differential cross-section w(0) = \A + 5A a + SAf,\ 2 one 
finds that Sw{9) = m[A*(5A a + SA b )}, which yields 



5r- x = 



ndcf 



3? / dpdp'dkdk'S E 5 



1 — nn' 



2n i v 2 "J " a " M 'k 2 - {k') 2 + ^0 

k.k' lm r- f- n , n W\ « « 1 ,fk,k' 



V(\k - p|)fcr^; - 2V(\p' - p\)kt ^ n F (k)n' F ( P ), (10) 



f,k.k' 



where fcr ppl = |Tr [s n us n '(k + S • k')us m ] , ^^p'p' = l^ 1 [s n usn'] Tr [(k + S • k.')us r 

I 



In Eq. (fTOj) . the integration is constrained by the 
momentum and energy conservation imposed by 8m — 
8 2 (p — p'+k' — k) and 8e — 8(p—p') and we use the nota- 
tions n = p/p, n' = p'/p, m = k/fc. Two form-factors of 
Bragg scattering, T and T represent the Fock and the Ha- 
tree contributions, respectively. Assuming x — y isotropy 
of disorder, we average the relaxation rate r _1 over the 
directions of the initial wave vector p. 

The Bragg scattering correction to the momentum re- 
laxation rate (|f Op can be separated into zero-temperature 
and temperature-dependent parts, <5t -1 = St^ 1 + 8t^ 1 . 
The linear temperature dependence of St^ 1 [3( is due to 
a singularity at k = k' in the integral (fT0|) . To evaluate 
its contribution we extend the analysis presented in Ref. 
3]. First, we use 5m and Se to rewrite the denominator 
in the integral (fTO]) as k 2 - [k') 2 = 2(k - p)(p - p') = 
2kp(cos9 - cos 9') - 2p 2 [l - cos(6» - 9% where 9 = kp 



and 9' = kp'. The integrand is singular at those points 
where the denominator vanishes, except for 9 = 9' where 
the singularity is cancelled by the factor [1 — cos(# — 9')]. 
For given k and p the locus of singular points is a contour 
in the plane (9, 9') which we show in Fig.[5]for various ra- 
tios k/p. As k changes from k < p to k > p, two parts of 
this contour coalesce creating a double pole at the point 
(9,9') — (0,7r). It is convenient to rewrite Eq. (fT0|) as 



5t^ = 



^def 
4 7r 4 w 2 



dpdkn F (k)n' F (p)g(k,p) — 5t q , 



where the function g(k,p) is the result of the integration 
of Eq. (fT0|) over p' and k' and the angular components of 



p and k. The leading term in St^ is linear in tempera- 
ture because of the jump, g(p + 0,p) — g(p — 0,p) = Ag(p) 
of the function g at k = p resulting from the double pole 
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a) 




b) 




FIG. 1: The figure is a schematic representation of the pro- 
cesses renormalizing the scattering amplitude in the Hartree 
a) and the Fock b) channel. 
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FIG. 2: The structure of singularities in the integration. The 
contours show the position of the singularities of the inte- 
grand in the 9, 9' plane. The arrow indicate the direction of 
increasing k/p. 



at (6,9') = (0,tt) for k = p: 
n dcf Ag(k F 



4tt 4 u 2 



T : 



-T^£^r y(o), (ii) 



u 2 s1 {1 + S sz ). 



- E 

s,l—x,y,z 



When deriving Ag(kp) in Eq. (|1 lj) . we took into ac- 
count that the contribution to the jump of g comes from 
backscattering, (9,9') = (0,7r), so that we substituted 
p = k and p' = p in all smooth angle-dependent 
functions. This produced the exchange term form- factor 
To = (Pp'-p) averaged over the incidence angle, and 
also resulted in the vanishing of the Hartree term due to 
the absence of backscattering of chiral electrons off the 
Coulomb potential, Tp'Ip = 0. 

Finally, we arrive at the following expression for the 
dominant temperature-dependent part of the resistivity, 



R(T) - R(0) 



h ATV(0) 
e 2 hv 2 epT t 



r . (12) 



It is interesting to compare the latter result with that 
derived |3[ for electrons in a simple-band 2D semicon- 
ductor or metal. In the latter case, R(T) is formed by 
the competition of the Hartree and Fock contributions 
which have different signs, which may lead to the change 
of the size and even the sign of the effect upon variation of 
the electron density or spin polarisation of the 2D gas by 
an external magnetic field. In graphene the T-dependent 
correction to the resistivity is only due to the exchange 
interaction, so that it is negative and its density depen- 
dence tracks the density dependence of the interaction 
constant V(0). For a 2D-screened Coulomb interaction, 
V(q) = 2n(e 2 /x)/(q + k), we estimate V(0) = ±7 _1 , 
which leads to the result in Eq. (fT]). The other difference 
between graphene and usual semiconductor structures is 
that in the former the linear in T correction is caused 
only by atomically sharp disorder (e.g., structural de- 
fects in graphene, chemical deposits, substitutional dis- 
order and contact with an incommensurate lattice of a 
substrate), whereas in the latter it is formed by all scat- 
terers. This means that the temperature dependence of 
resistivity should be weak (or even absent) in a suspended 
graphene sheet or a sheet loosely attached to the sub- 
strate, and that it can be pronounced in devices where 
graphene is strongly coupled to the substrate. 

In conclusion, we described the Friedel oscillations in- 
duced by scatterers in graphene and their effect on resis- 
tivity. The general form of FO is given in Eq. |8]). One of 
the immediate implications of this result is a r~ 2 depen- 
dence of the RKKY interaction between two magnetic 
impurities in graphene. This is because a substitution 
atom with spin S will generally create a perturbation 
(S • s e )u for the electron spin s e with u containing all 
symmetry-breaking terms in Eq. j3|. We also showed 
that due to the chirality of graphene electrons, the lin- 
ear temperature correction to the resistivity [ Eqs. flTJ) 
and (TP2|) ] is caused by atomically sharp defects rather 
than by Coulomb charges in the insulating substrate, so 
that its measurement can be used as a tool to test the 
microscopic composition of disorder. 

We thank I.Aleiner and A.Ludwig for discussions. This 
project was funded by the EPSRC grant EP/C511743. 



[1] J. Friedel, Phil. Mag. 43, 153 (1952); 
[2] K.H. Lau and W. Kohn, Surf. Sci. 75, 69 (1978) 
[3] G. Zala, B.N. Narozhny, and I.L. Aleiner, Phys. Rev. B 
64, 214204 (2001); ibid. 64, 201201 (2001); ibid. B 65, 
020201 (2002) 

S. Das Sarma and E. H. Hwang, Phys. Rev. Lett. 83, 164 
(1999); Phys. Rev. B 69, 195305 (2004) 



5 



[5] 
[6] 
[7] 

[8] 
[9] 
[10] 



[11] 
[12] 



[13] 
[14] 
[15] 



A.M. Rudin, I.L. Aleiner, and L.I. Glazman, Phys. Rev. 
B 55, 9322 (1997) 

F. Stern, Phys. Rev. Lett. 44, 1469 (1980); A. Gold and 

V.T. Dolgopolov, Phys. Rev. B 33, 1076 (1986) 

Y.Y. Proskuryakov et al, Phys. Rev. Lett. 89, 076406 

(2002); Z.D. Kvon et al, Phys. Rev. B 65, 161304 (2002); 

A. A. Shashkin et al, Phys. Rev. B 66, 073303 (2002); 

V.M. Pudalov et al, Phys. Rev. Lett. 91, 126403 (2003); 

S.A. Vitkalov et al, Phys. Rev. B 67, 113310 (2003) 

K.S. Novoselov et al., Science 306, 666 (2004); 

K.S. Novoselov et al., Nature 438, 197 (2005) 

Y. Zhang et al, Nature 438, 201 (2005); Y. Zhang et al, 

Phys. Rev. Lett. 94, 176803 (2005) 

R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical 
Properties of Carbon Nanotubes, Imperial College Press, 
London 1998 

T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005) 
K. Nomura and A.H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006); T. Ando, J. Phys. Soc. Jpn. 75, 074716 
(2006) 

E. McCann, et al, 



£(|7ra 1 , 0), where £ = ±1 and a is the lattice constant. 
In the basis [</>k + ,a, <)>k , 



t>K._ B, <pK._,A\ 



time rever- 



cond-mat/0604015 
■B. Efetov, |cond-mat/0607200| 



I.L. Aleiner and K 

Corners of the hexagonal Brilloin zone are 



sal, T(W) of an operator W is described by T(W) = 

(ILj ® o~x) W*(H X <S> Ox). This can be used to show that 

T(£ s ) = -E s , T(A,) = -A,, and T(E S A ; ) = E a A ; . 
[16] E. Fradkin, Phys. Rev. B 33, 3257 (1986); E. McCann, 

V. Fal'ko, Phys. Rev. B 71, 085415 (2005); M. Foster, A. 

Ludwig, Phys. Rev. B 73, 155104 (2006) 
[17] S.V. Morozov et al, Phys. Rev. Lett. 97, 016801 (2006) 
[18] T. Ando, T. Nakanishi, and R. Saito, J. Phys. Soc. Japan 

67, 2857 (1998) 
[19] V. Cheianov and V.I. Fal'ko, Phys. Rev. B 74, 041403 

(2006) 

[20] T. Ando, A. Fowler, F. Stern, Rev. Mod. Phys. 54, 437 
(1982) 

[21] We neglect the non-conservation of the sublattice and 
valley indices in the e-e scattering. The corrections due to 
such processes are small, as V(2n /a) /V(kp) ~ k?a <C 1, 
where a is the lattice constant. 

[22] The 1/r 2 FO in a gas of relativistic 2D fermions vanish 
in the massless limit. See: D.H. Lin, Phys. Rev. A 73, 
044701 (2006) 



